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Abstract 

We derive expressions for the form factors of the quantum transfer matrix 
of the spin- 5 XXZ chain which are suitable for taking the infinite Trotter 
number limit. These form factors determine the finitely many amplitudes 
in the leading asymptotics of the finite-temperature correlation functions 
of the model. We consider form-factor expansions of the longitudinal 
and transversal two-point functions. Remarkably, the formulae for the 
amplitudes are in both cases of the same form. We also explain how to 
adapt our formulae to the description of ground state correlation functions 
of the finite chain. The usefulness of our novel formulae is demonstrated 
by working out expHcit results in the high- and low-temperature Hmits. We 
obtain, in particular, the large-distance asymptotics of the longitudinal two- 
point functions for small temperatures by summing up the asymptotically 
most relevant terms in the form factor expansion of a generating function 
of the longitudinal correlation functions. As expected the leading term in 
the expansion of the corresponding two-point functions is in accordance 
with conformal field theory predictions. Here it is obtained for the first time 
by a direct calculation. 
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1 Introduction 

For a long time form factors have been the main tool for studying the correlation 
functions of integrable massive quantum field theories through their Lehmann represen- 
tation [20,47]. For these theories they can be determined as solutions of appropriate 
functional equations [48,49]. 

Form factors of integrable non-relativistic models can be studied by means of the 
vertex operator approach [18] and by means of the algebraic Bethe ansatz [23-25, 27, 
28, 46]. In this work we apply the latter approach in order to derive expressions for 
the form factors of the quantum transfer matrix [51,52] that are suitable for taking the 
limit of an infinite Trotter number. Similar formulae hold for the form factors of the 
usual row-to-row transfer matrix for finite chains. In both cases the form factors will 
be described in terms of solutions of linear and non-linear integral equations, similar 
to those appearing in the description of the thermodynamics [29, 30] and correlation 
functions [4, 14] of the model. We shall call such form factors thermal form factors. 

At finite temperature the correlation functions of the XXZ chain decay exponentially. 
Their asymptotic behaviour is described by correlation lengths, determined by the 
ratio of one of the sub-dominant eigenvalues of the quantum transfer matrix and 
its dominant eigenvalue, and by amplitudes given as products of two corresponding 
thermal form factors. While the correlations lengths were studied thoroughly in the 
literature [29, 30, 33,43], so far rather little attention has been paid to the amplitudes. 

Our expressions for the form factors involve, except for integrals over functions 
obtained as solutions of integral equations, also Fredholm determinants of the cor- 
responding integral operators. We consider form factors arising in expansions of a 
generating function of the longitudinal correlation functions and of the transversal 
two-point correlation functions. The usefulness of our formulae is demonstrated in 
several limiting cases for high and low temperatures. 

The analysis of the low-temperature hmit is delicate and cumbersome. In this work 
we concentrate on the longitudinal correlation functions in the critical regime. The 
low-temperature analysis of the transversal correlation functions in the critical regime 
is similar but requires a separate calculation which we leave for future work. Quite 
generally, critical models at low temperature are characterized by the fact that infinitely 
many of their correlation lengths diverge [10, 11]. We shall see that the corresponding 
amplitudes vanish algebraically as functions of the temperature with the same critical 
exponents that govern the spatial decay of the ground state correlation functions. The 
situation is very much like for the finite-size dependence of the 'critical amplitudes' 
occurring in the form factor expansion of the ground state correlation functions [23, 25]. 
A similar temperature dependence was also reported for the amplitudes in a series 
expansion of the density-density correlation functions of the Bose gas with delta- 
function interactions [37, 38] obtained from a multiple-integral representation of the 
generating function. 

In fact, we shall closely follow the work [37] in our low-temperature analysis. This 
allows us to keep the presentation short, as many details can be proved along the lines 
of [37]. As in [37] we can also sum up the critical particle-hole contributions to the form 
factor expansion by means of a summation formula proved in [2 1 , 24, 40] . As a result 
we obtain the low-temperature large-distance asymptotics of the longitudinal correlation 
functions. The form of the exponential decay is in accordance with conformal field 
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theory predictions [10, 1 1]. The ampUtude of the leading oscillating term so far could 
not be obtained by approximate field theoretical calculations. It is the same as the 
ground state's leading oscillatory amphtude which was obtained by a direct asymptotic 
analysis of multiple integrals in [22]. 



2 Hamiltonian and correlation functions 



The XXZ-chain in a longitudinal magnetic field h is defined by the Hamiltonian 



j=-L+l / Z , 



(1) 



j=-L+l 



acting on 



. Here the aj, a = x,y,z, are the Pauli matrices and periodic boundary 
conditions, a"^ = a", are implied. J and A are the exchange energy and the anisotropy 
parameter. We shall employ the standard parameterization A=(^ + ^^)/2,^ = e^. 
The Hamiltonian (1) is integrable with /?-matrix 



R{X,/j) 



b{X,/j) 



/I 0\ 

b{X,n) c{X,iS) 

c{X,ii) b{X,n) 

Vo 1/ 



sh(X.— /i) 

SKX-Zl + Tl) ' 



c{X,n) 



sh(Ti) 



sh(X-/i + Ti) 



(2a) 



(2b) 



For the exact calculation of thermal averages [14] at temperature T we have to 
express the statistical operator of the grand-canonical ensemble 



(3) 



in terms of the /^-matrix (2a). For this purpose we introduce an auxiUary vertex model 
with monodromy matrix 



P 



(4) 



It acts non-trivially in the tensor product of an 'auxiliary space' j = —L + 1 , . . . ,L and 
N elN 'quantum spaces' l,...,N which are all isomorphic to C^. The superscript ti 
denotes the transposition with respect to the first of the two spaces R(k,ij) is acting on. 
The parameters P and k are related to the temperature T and to the magnetic field h, 



P = 



2ysh(Tl) 



K = 



2t[\T 



(5) 



Tj{X) is a 2 X 2 matrix in auxihary space j. Its matrix elements satisfy the Yang-Baxter 
algebra relations 

Rjk{X,ix)Tj{X)Tk{jS) = TkiM)Tj{X)RjkM (6) 



by construction. 
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Using the monodromy matrix (4) we can express the statistical operator (3) as a 
Umit, 

PL = lim Tri. jv{r-z.+i (0) . . . 7i(0)} . (7) 

We call this limit the Trotter limit and the number the Trotter number. For finite 
Trotter number the operator 

pN,L = TT-,.M{T.L+liO)...TLm (8) 

approximates the statistical operator. It can be used to define finite Trotter number 
approximants of correlation functions of m local operators acting on sites 

1 , . . . , m of the infinite chain 



<0(/) . . . = hm TrWp...Or--0;n^} 
^ '^^ L^°° Tr_z.+i...z.{Piv,L} 



Trj jv{Tr^{r(0)}Tr{O(i)r(0)} . . .Tr{OWr(0)}Tr^-'"{r(0)}} 
Tri...]v{Tr2^{r(0)}} ' 

Here the expression 

Z;v,L = Tri...^{Tr2^{r(0)}} (10) 

in the denominator is the finite Trotter number approximant to the partition function 
Zl = Tr ^,+1 x{exp(— ///r)} of the XXZ chain of length 2L. The operator occurring 
on the right hand side of equation (10), 

t{X)=TT{T{X)}, (11) 

is called the quantum transfer matrix [51]. It has the pecuMar property that the spectrum 

of t{0) contains a unique real eigenvalue Ao(0), whose modulus is strictly larger than 
those of all other eigenvalues, even if the Hamiltonian (1) is in the critical regime | A| < 1. 
This implies that t{0)^ in (9) acts as a projector onto the corresponding eigenstate I^I'o)^ 
when L becomes large. It follows that 

(.H\ _ {"Vol Tr{0(i)r (0)} . ■ .TrjoWr (0)}|>yo) 
^ ' /AT- (*o|*o)A^(0) • ^'^^ 

We shall call |*Fo) the dominant eigenstate and Ao(A,) the dominant eigenvalue. We see 
that the dominant eigenstate determines all thermal correlation functions of the system 
in the thermodynamic limit, or, to phrase it differently, that it determines the state of 
thermal equilibrium completely. The physical correlation functions follow from (12) in 
the Trotter limit N ^°o. For more details see [14]. 



3 Bethe ansatz and auxiliary functions 

Eigenvectors and eigenvalues of the quantum transfer matrix can be constructed by 
means of the algebraic Bethe ansatz. This follows since the monodromy matrix elements, 
for which we use the standard notation 
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satisfy the Yang-Baxter algebra relation (6) and since a pseudo vacuum |0) = [( j ) (8) 
(o)] ^'^^^ exists on which the monodromy matrix elements act as 

C{X)\0)=0, 

A{Xm=q^b^-l^,X)\0)=a{X)\0), (14) 

D{xm=q--b^x4)\o)=d{xm. 

Then (see e.g. [36]) 

\Wn)=B{lM).-.B{h)\0) (15) 

is an eigenvector of the quantum transfer matrix t{X) if the Bethe roots Xj, j = 1,... ,M, 
satisfy the system 

of Bethe ansatz equations. The corresponding eigenvalue is 

A„{X) - a{X) [[ +dix) [[ ,h(x-x,) • ^^^^ 

A left eigenvector of the quantum transfer matrix with the same eigenvalue is 

{^n\ = {0\C(ki)...C{XM), (18) 

where {0| = |0)' is the left pseudo vacuum satisfying {0\B{7i) = 0. 

We have used « = 0, 1, ... ,2^ — 1 to enumerate the eigenstates. This means we 
assume implicitly that, for generic q and non-zero K, the sets of right and left eigenvec- 
tors defined by (15) and (18) form bases of (C^)*^^ and of its dual. We would like to 
emphasize that this point is important for the interpretation of our results below, but 
not so much for the results themselves. The quantum transfer matrix of the XXZ chain 
can be seen as a special inhomogeneous case of the usual transfer matrix of the model 
(see e.g. [4]). For a discussion of the issue of completeness of the Bethe states of the 
latter, as far as it is relevant in the present context, the reader is referred to [27]. In the 
following we shall assume that the eigenvalues are ordered by size in such a way that 
Ao(0)>|Ai(0)|>|A2(0)|>.... 

Along with the quantum transfer matrix t{X) we consider a 'twisted version' t{X\a) 
obtained from t(k) by shifting K -> K-|-a. Then a{X) — )■ q^a{X), d(k) — )■ q~^d(k) and 
the Bethe roots are altered as well. For the corresponding eigenvalues and eigenvectors 
we use the notation A„(k\a) and \^'^). 

Our goal is to write the form factors of the quantum transfer matrix in a form that 
admits the Trotter limit. For this purpose we introduce certain auxiliary functions 
generalizing some of those functions that proved to be useful in the description of the 
finite-temperature correlation functions of the model [4, 14]. For every regular solution^ 
of the Bethe equations (16) we define an auxiliary function 



^We call a solution of the Bethe equations regular if it corresponds to an eigenvector of the form (15). 
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In the twisted case we write o„(A,|a). We also define 

P.(.|a).--^. 

In the Trotter limit the functions a„(k) are uniquely determined by the integral 
equation 

ln(a„(?i)) = -(2K+Af-2M)ii-pe(X)- / ^K{X-ij)ln(\ + an(jj)) . (21) 

Je„ 2711 

Here we have introduced the kernel 

K(k)=Ko{X), Ka(X)=q-''cth{X,~^)-q''cth{X + ^), (22) 
which is the derivative of the bare two-particle scattering phase, and the bare energy 

e(k) = cth(X) - cth(^ + 'n) . (23) 

In equation (21) the information about the specific excitation is hidden in the choice 
of the integration contour C„. For the dominant state all Bethe roots are located inside 
the strip between —v\/2 and v{/2!^ All zeros of Ao(?i|a), which are zeros of 1 + ao{X\a) 
as well, are located outside this strip. These zeros are called holes. Thus, the contour Co 
of the dominant state consists of two straight Mnes parallel to the real axis at ±ti/2 and 
encircling the real axis in positive direction. If X is on the contour we infinitesimally 
deform it in such a way that the point at the opposite side of the contour (k+r{orX — r\) 
is outside. Excitations are characterized by a finite number of holes X'j inside Co and a 
finite number of Bethe roots X^ outside. The contours C„ are obtained by deforming So 
in such a way that the X^j are excluded from C„ but the X^ are included (see figure 1). 
For our convenience we shall also assume that all Bethe roots of the dominant state are 
located inside S„ while the corresponding holes remain outside." 

The contours Q„ can be straightened. Then modified driving terms depending 
explicitly on the particle and hole parameters X^ and X'j appear, which can be determined 
by a set of subsidiary conditions of the form a„(X[^) = — 1, a„{X'j) = — 1 (see [30]). Such 
type of formulation is useful for numerical calculations and for the low-temperature 
analysis of section 8. For the derivation of form-factor formulae, however, the deformed 
contour formulation seems to be more convenient. 

The functions p„(X|a) can be represented as integrals over the auxiliary functions, 

p.(,„,.,"^«e.p{/Je„-.),„(l±^)}. (M, 

where X is located inside the contour S„. The number A^/2 — M is the eigenvalue of the 
conserved z-component of the pseudo spin, 11" = ^ I-^=i (~^)^^)- '^^^ discussion of 
the form factors below we have to distinguish two cases: (i) A'^/2 — M = 0, which is 
relevant for the longitudinal correlation functions, and (ii) N/2 — M=l, which applies 
to the transversal correlation functions. In both cases the right hand side of (24) is well 
defined and non-trivial in the Trotter limit. 

""For definiteness we describe the situation in the critical case < A < 1 corresponding to t] = — iy, 
Ye(0,3t/2). 

II This implies that we exclude the free Fermion case j=n/2 from our considerations. 
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Cn ^ 

-iY/2 

Figure 1: Sketch of a contour C„. Hole parameters inside the 'physical strip' between 
—iY/2 and iY/2 are excluded from the contour, and particle parameters (Bethe roots) X.^ 
outside the physical strip are included into the contour. 



4 Form factor expansion of correlation functions 

We shall consider form factor expansions of the longitudinal and transversal two-point 
functions (cJia^_,_i) and For the longitudinal two-point functions we start 

with a generating function [17] which is closely related to the twisted transfer matrix. 
Setting 

1 m 

5H = -£a} (25) 

the generating function is defined as ^^^"^('")^, where a is a complex parameter. In 
fact, 

(^^^^+l) = 2^^35a<^'"'('"+^)>[^„, (26) 

where is the difference operator defined by Dmf{m) = f{m) — f{m — 1). Note that 
the generating function provides more information. The one-point functions are given 
by 

(a^)=D^^^„(^2oc5(m+l)^^ ^ (27) 

a=0 

and the emptiness formation probability can be obtained as the limit 

{eiX...e^X)= hm ^-'"«(^2<x5(m)^> _ (28) 
^ ' ReTia->+<>° ^ ' 



Here ej^ is a projector onto the local spin-up state. 
It follows from (12) that 
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is the finite Trotter number approximant to the generating function. Expanding 
?'"(0|a)|*I'o) in the basis of Bethe states {I*!*")} we obtain the form factor expansion 

{d'^^^-X = I A„(a)p:(0|a) , A„(a) = ■ (30) 

Here we have adapted the numbering of the states. We count only those Nm = (^), 
M = N/2, states which have non-zero overlap with the dominant state. 

Inserting now (30) into (26) we obtain the form factor expansion of the connected 
longitudinal correlation functions, 



Nm-I , 1 _i\2 

(a^ic4+i)iv- (al)iv(a^+i)Ar = I A- [pl-pn') . (31) 



Nm-I 

I 

n=l 

In this expression we introduced the abbreviations 



P, = P.(0|0). ^J = l3;;.^.(„)|^„=|^{||i|M. <32, 

By and we mean the derivatives of the corresponding a-dependent vectors 
with respect to Tja at a = 0. 

Note that 1 > |pi | > jpzl > • • • by construction. This means that (31) is the finite- 
temperature asymptotic expansion of the longitudinal correlation functions. The leading 
asymptotics is determined by the first few terms. The correlation length of the longitu- 
dinal correlation functions is 

^i = -Ari- (33) 
ln|pi| 

This correlation length was studied in [13,30,33,43,53]. So far the amplitudes in (31) 
were only studied nimierically [13] for finite Trotter numbers. Below in section 5 we 
shall calculate them analytically in the Trotter limit. 

Using once more (12) we obtain a finite-temperature asymptotic expansion for the 
ttansversal correlation functions, 

K<+i)/v=I:a„+p™, (34) 

n=l 

where 

+ ^ (^o|g(0)|*l^„) (%|C(0)|%) 
" A„(0)(»Fo|»Fo) Ao(0)(%|^„)- ^ ^ 

Again the sum is restricted to the non-vanishing amplitudes (for which necessarily 
M = N /2 — \) and the numbering is adapted (in this case we start with '1', since we 
have reserved the label '0' for the dominant eigenstate). Accordingly the p„ in the above 
formula (34) belong to excitations with M = N /2 — 1 Bethe roots and are different 
from the p„ in equation (31) belonging to excitations with M = N /I Bethe roots. As 
before the correlation length is defined by (33) with the appropriate pi inserted. The 
amplitudes in the Trotter limit were heretofore unknown and will be calculated in 
section 6. 
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5 Amplitudes in the asymptotic expansion of the longitudi- 
nal correlation functions 

The basic equation underlying the algebraic Bethe ansatz approach to the correlation 
functions of the XXZ chain is the following scalar product formula due to Nikita 
Slavnov [46]: 

(0|C(/ii) . ..C{mm)B{Xm) . ..B{Xi)\0) = (0|C(Xi) . . .€^8^ • • .B{tii)\0) 



M M 1 



detM(e(X,y -;Ujt) - e(^/: - X^-)o„(/ifc)) 



det, 



(36) 



Here the Xj, j = 1, . . . ,M, are a solution of the Bethe equations (16) and o„ is the 
associated auxiliary function. The /uj in this formula are free. They may or may not be 
a solution of the Bethe equations. In particular, we can take the limit /^j — > Xj in (36) 
and obtain the 'norm formula' for Bethe states [35], 



M 



1-7=1 



det<^ di - 



K{Xj — Xk) 



(37) 



Using (36) and (37) in (32) or (35) we obtain expressions for the amplitudes in 
terms of Bethe roots (compare [28]). These are not directly suitable for taking the 
Trotter Umit. Still, after a sequence of elementary manipulations (see appendix A) we 
obtain the following result. 



A„(a) 



^ p„(X,|a) l dew{8rg^3C_«(^,-^.)} 
,yP.(A^7icc)J detM{5^^XCu,-^,)} 

detM{5^-^3C(X,--X,)} 



(38) 



In this equation M = N/2, the Xj, j = l,...,M, are the Bethe roots of the dominant 

state, while the i^j, j = 1, . . . ,M, are the Bethe roots of an excited state of the twisted 
transfer matrix t{X\a). The kernel functions in the determinants are defined as 



X{X)=Xo{X), Xa{X) 



g(a-l)(X-Ti) Q(a-l){X+r\) 

sh(x.-'n) sh(x.+'n) 



(39) 



The right hand side of (38) is suitable for taking the Trotter limit. For the product 
this is obvious from the representation 



^ Pn{Xj\a) 

1 1 / I X = exp 



dX p'„{X\a) /\ + an{X\a] 



e„27ii Pn{X\a) V l + ao(>^) 



(40) 



The determinants are all of the same structure. Expanding e.g. the first determinant in 
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the numerator we obtain 



_^ y P«(i";l«) P«(A'/fc|«) 

= 1 + ^^ 

<:=1 



X-a(O) 
3C_a(/^^-)Uj) X_a(0) 



n/" dm«(v,)jdet{3C_a(v,-v^)}, (41) 



where we have introduced the 'measure' 



dA,p„(A,|a) 



2jti(H-o„(?i|a)) ■ 



(42) 



In the Trotter limit A'^ = 2M — oo the right hand side of (41) converges to the Fredholm 
determinant of the integral operator %-a defined by the kernel %-a, the measure dm" 
and the contour C„, 



lim detls^ - P;iM^3C_a(;U,- -/i,)) 



= 1 + 1 



n/ dm«(v,)ldet{DC_a(v,-v«)}= det {1-0C_„}. (43) 



The other determinants in (38) can be treated in a similar way. Introducing the 
measures 



dm(X.) = 
dmU^) = 
dm«(X) 



dX 



2jii(l + oo(X)) ' 
dX 

27ti(l + o„(X|a)) ' 
dXp-\X\a) 



2ra(l + oo(^)) 
we obtain 

^"^"^"'''n427rip„(X|a)^l 1 + aoW J i 

deta„,» ,e„ { 1 - X-g} deta,„.x .e„ { 1 - Xg} 

X ^ ^ 

detd^a e„ {^-^} detdm,e„ {^-^} 



(44a) 
(44b) 
(44c) 



(45) 



for the amplitudes in the Trotter limit. At least at first sight this looks rather appealing. 
The amplitudes are entirely described in terms of functions which appeared earlier 
in the description of the thermodynamic properties, the correlation lengths and the 
correlation functions of the model. These are the auxiliary functions a„ (see [30]), the 
eigenvalue ratios p„ (see [4, 8]) and the deformed kernel %a. 

The latter is a close relative of the deformed kernel K^,, equation (22), which first 
appeared in the context of the description of the correlation functions of the XXZ chain 
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in [4, 7]. It seems that the kernel Xa is more suitable in the analysis of the scaling limit 
towards conformal field theory [2, 8], while Ka seems to be more convenient for the 
high-temperature analysis. This is also reflected in the different symmetry properties of 
the two kernels, 

aCa(-X) = X2-aW, Ka{-X)=K_a{^). (46) 

Since the two kernels are related by 

e-«^3C„(X) = Ka{X) - + ^« , (47) 

the integral operators in the determinants in (45) differ from the corresponding integral 
operators with kernel only by a one-dimensional projector. Hence, using also the 
fact that %o = Ko = K, we can rewrite (45) as 



A„(a)=a;a:exp ^^^In 



dX p',Qi\a) r \ + an{X\a) 
2m Pn(X\a) °V l + cio(5^) 



(48) 



^ detdM« ,e„ { 1 - ^-a} detd^a { 1 - ^a} 
detd^a e„ { 1 - ^} detdm,e„ { 1 - ^} 
where 

o^= hm c+{X), GZ= hm a_(?i) (49) 

ReX.— >— oo ReA,— oo 

and where the two functions a+(k), 0-(X) are generalizations to excited states of the 
generalized dressed charge functions introduced in [5]. They are defined by the integral 
equations 

a+(X) = 1 + / dmliM)Ka{X -/i)a+(;u) , (50a) 
G-{X) = \+ f dm" (Ai)a_ (jii)Ka i/n-X). (50b) 

Using the techniques developed in [5] it is not difficult to show that 

_2 f /■ dX f l + a„(X|a 



alal = ch-"<^ / — In ' ' } . (51) 

When both sets of rapidities in (36) satisfy the Bethe equation we obtain two 
different formulae by using either of the two corresponding auxiliary functions on the 
right hand side of the equation. In order to derive (45) we used both possibilities in a 
symmetric fashion and ended up with a particularly symmetric result. Using two times 
the same formula and also (37) we obtain alternative expressions for A„(a) containing 
the square of one of the two determinants that appear in the numerator on the right hand 
side of (45) and a different prefactor. Comparing these asymmetric formulae with (45) 
we obtain the following interesting relation between the Fredholm determinants in the 
numerator, 

detd^«,e„{l-Xa} _(M^_f dX^^fl + aniXla) 



detd^a_e„{l-3<:-oc} I Pn(0|a) Je„ni \ l + ooiX) 
dX 



e„27ti 



^l„(l+^W)-^Kl+a,,.|a)) 



(52) 
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Equations (45) and (48) are the main results of this section. They are valid for 
finite Trotter number as well as in the Trotter limit. For finite Trotter number the right 
hand sides of these equations are represented by finitely many integrals (see (41)). We 
will show below that the high- and the low-temperature limits of the amplitudes can 
be extracted from these formulae. Their suitability for numerical calculations will be 
explored in future research. 



6 Amplitudes in the asymptotic expansion of the transversal 
correlation functions 

In this section we consider a slight generaUzation of the amplitudes A^+, equation (35), 
depending on a spectral parameter ^ and on a twist parameter a. We define 

and 

A-+(^)=F_(^)f+(^). (54) 
For F+(^) we obtain in appendix B 

G;(^)e-^nf^i^(^y)e^^ detM-i{6i-g|;|gg/:i_.(^,-^.)} ^^^^ 

Here 

G;(^)= hm G+(X,^), (56) 
and G+ {X, ^) is the solution of the linear integral equation 

G+(X,^) = -cth(X-^) 

+ ^«-Vn(^|a)cth(^-^-ri)+ / dm«(A/)i^i_a(X-A')G+(A',^)- (57) 

The simplest way to obtain a formula for the amplitudes from (55) is to combine 
(55) with (36) and (37). Then 



G;(^) ' ^.-I-.V.-^) nfj/ AolA-i) 



Pn(^|a)nf=/^-«^(/^,)PnC",|a) nf=iAo(:^,-) 

detM-i{6^ - -A^^)} detM{8i - ^^K{Xj-Xk)] 



(58) 
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As before the determinants turn into Fredholm determinants in the Trotter limit. Some 
more care is necessary when dealing with the prefactors. All in all we obtain 

dX p'„{X\a) ( l + a„(X|a 



where 



l + Oo(X) 

7 (59) 

detd;„a e„ 1 1 - ^1 detd«,e„ i^-J^) 



^1 p„(0|a) ye„ 711 ^ V 1 



+ 



/e.li[^'"(..MX))-^K...(X|a)) 



Oo(?l) 

(60) 



Note the similarity between the exponential term on the right hand side of this 
equation and the right hand side of equation (52). This similarity hints that a formula 
hke (52) might also exist for the present case and that a simpler, more symmetric 
formula for A~+(^) can be derived. 

There is indeed an alternative, less obvious formula for the form factor F_ (^), 



(61) 



where 



Gl{k)= lim G-iXX). (62) 
ReA— 

and G_ (X., ^) is the solution of the linear integral equation 

G_(X,^) = -cth(X-^) 

+ ^«+ip„-i(^|a)cth(X-^-il)+ / dm«(;u)G_(Ai,^)i^i+aC"-:^). (63) 

The derivation is shown in appendix B. 

Combining (55) and (61) we obtain a more symmetric expression for the ampMtudes 

n (^i+a_^-l-a)(^a_^-a)p^M-lp^(^.|„) 

detM- 1 1 - " A'^) | detM [^^ " | 
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In the Trotter limit it turns into 



,a+K-l I -a-K+1 




detd^a ,e„ { 1 - ^i-a} detd;„a e„ { 1 - Ki+a] 
detd^a {\-K] detd;„,e„ {^-K] 



(65) 



What we find remarkable about this expression is that it is basically of the same 
form as the expression (48) for the amplitudes in the form factor expansion of the 
generating function of the zz-correlation functions. The exponential term and the ratio 
of the products of two determinants appear in both cases, and the prefactor is a product 
of two functions satisfying 'conjugate' linear integral equations. All expressions are 
parameterized by the functions p„ and a„ which depend upon the specific excitation 
under consideration and are, of course, different in both cases. In appendix B we argue 
that limo(_j,o (^) = 0, such that the prefactor in (65) is non-singular in this limit. 

Following [8, 9] we can introduce 'the spin' of a local operator as the eigenvalue of 
the z-component of the operator of the total spin when acting adjointly on it. Then (5 , 
a^, a+ have eigenvalues 5 = — ,0,+, and we can interpret the parameterization of the 
kernels in the Fredholm determinants in the numerators of (48) and (65) as ^T^-a and 
Ka-s- In this combination the spin also appears in the TQ equation. 

7 Generating function at high temperatures 

We are convinced that our formula for the ampHtudes are not only aesthetically appealing 
but that they will also prove to be efficient for the future analytical and numerical 
analysis of the long-distance asymptotic behaviour of the finite-temperature correlation 
functions of the XXZ chain. In this and in the following section, we shall support this 
claim with two examples concerning the high- and the low-temperature asymptotics of 
the generating function of the longitudinal correlation function. 

First of all, we use (48) to calculate the leading high-temperature contribution to 
the generating function {^q^'^^^'"^^ Since the next-to leading eigenvalues vanish in the 
high-temperature limit, we only have to consider the first term n = of the expansion 



As is well known (see e.g. [55]), the dominant eigenvalue of the (twisted) quantum 
transfer matrix and the corresponding auxiliary function can be expanded in a power 
series in l/T. Let us briefly review this procedure. 
It follows from (21) with M = N /2 that 



if P = K = 0. In order to calculate the next order, we assume the asymptotic expansion 



(30), 



{q^'^^"''^) =Ao(a)pS'(0|a) + 0(l/r'"). 



(66) 



ao(X|a) = q 



(67) 




(68) 
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By substituting (68) into (21) and expanding in powers of l/T, we see that ai (k\a) is 
the solution of the linear integral equation 



(69) 



ai im = -h- 27 sh(ii) e(X) - J^^ ^.K{X - ,u) . 

Obviously, the function 

a,(X\a) = -/i-27sh(ii) (^e(X) - (70) 

is a solution of (69). Since ao(^) = Oo(X|0), the corresponding expression for oo(X) 
can be obtained from (70) by setting a = 0. These results can be substituted into (24), 
after evaluating the contour integrals we arrive at 

PoiMa) = ^ ^ 1 + — — + — — K(X) +0(\ T^). 

' ^ 2 I 2T q^ + q-^ IT \q«- + q-'^ J ^ 'j^ ^ ' > 

(71) 

By the same reasoning, we obtain 

a;a- = l + 0(l/r^), expj/ ^ ^^^f 1±^H = 1 + o(i/r^) . 
+ " V/ ^\yeo27tipo(X|a) Vl + ooWyj ^' ' 

(72) 

Lastly, we have to consider the Fredholm determinants appearing in (48). Since every 
pole of the measure comes with a factor l/T, we observe that the A:-th term of the series 
(43) is of the order 0(1/7^^). Hence, for our purpose, it is sufficient to calculate the 
first term of the expansion (43). Using (70) and (71) we obtain after some calculations 

Ao(a) = 1 + ^ f ^^'i ' + 0(l/r2) . (73) 



T yq'^ + q-^ 
Combining (71) and (73) with (66), we end up with 

.f„.,/4(tz±lV(t^y-\0iVT^). (74, 



Note that this result was akeady obtained in [45] within the multiple-integral approach 
to correlation functions. 

Using (26) we can calculate the longitudinal correlators, 

(^i^2) = -Y + 0(l/7^'), (75a) 



(a^a^+i) = 0(l/r^), m>2. (75b) 
The emptiness formation probabiUty is given by 

(eil . ..e.l) = (0 (l + - I'" - l)y ) + 0(1/7^^) • (76) 

Note that (74), (75) and (76) are valid for all A > — 1 . It will be interesting to extend 
the high-temperature expansion systematically to higher orders. 
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8 Generating function at low temperatures 

The low-temperature limits of the massless and the massive regimes (|A| < 1 and A > 1, 
respectively) have to be analyzed separately. For A > 1 the spectrum of the quantum 
transfer matrix remains gapped as T — > 0, but for | A| < 1 it becomes gapless, and the 
long-distance asymptotics of the static correlation functions is expected to be described 
by conformal field theory. Here we consider the massless case at finite magnetic field, 
h> 0. The low-energy (large-distance) physics of our spin chain has then much in 
common with the physics of free Fermions. It is largely determined by particle-hole like 
excitations around the Fermi edge. As we shall see, the low-temperature analysis can 
be based on a generalization of the Sommerfeld expansion [50] to the interacting case. 
Such an analysis was partially carried out in [34], but here we have to go further and 
include the quantum numbers of the excited states. In the following we shall set r\ = — iy 
and, for simpUcity, impose the restriction < y < Jt/2. Then < A = cos(y) < 1. 

Since the spectrum becomes gapless for T — > 0, infinitely many terms in the form 
factor expansion contribute to the large-distance asymptotics of the correlation functions. 
As we shall see, each individual contribution vanishes as a fractional power of the 
temperature. The situation is very similar as in case of the so-called critical form factors 
defined between eigenstates of the ordinary transfer matrix [25] which vanish in the 
thermodynamic limit as negative fractional powers of the system size. We shall also 
observe a close analogy of our result with that obtained for the Bose gas by one of 
the authors [37, 38] using a different method. As in [37] we will utilize a remarkable 
summation formula, obtained in [21,24,40], which allows us to sum up the leading 
low-temperature terms. This is necessary since the limit involved in the infinite series 
of the leading terms and the limit T — ^ do not commute. 

The technique we are going to apply can be used for analyzing both, the longitudinal 
as well as the transversal correlation functions at small temperatures. Here we shall 
restrict ourselves to the analysis of the generating function of the longitudinal correlation 
functions and postpone the transversal case to a future publication. We start from the 
form factor expansion (30). In order to determine the leading low-temperature behaviour 
of the amplitudes and eigenvalue ratios we use (45) and (24). The full calculation is 
rather lengthy. It consists of the following steps. 

(i) Straightening of the contours in the non-hnear integral equations (21) and in the 
expressions for the amplitudes (45) and eigenvalue ratios (24). 

(ii) Low-temperature analysis of the non-linear integral equations. 

(iii) Low-temperature analysis of the eigenvalue ratios. 

(iv) Low-temperature analysis of the amplitudes. 

(v) Summation of the leading terms. 

Fortunately, the most sophisticated step (v) can be adopted from [24], and essentially 
only steps (i)-(iv) remain to be done. 



17 



8.1 Straightening the contours 

The contour Cq consists of two straight lines parallel to the real axis and intersecting 
the imaginary axis at ±Ti/2 such that the real axis is 'encircled' in counterclockwise 
direction. The contours C„ are deformations of Co such that all Bethe roots are included 
and all holes are excluded. To straighten C„ means to deform it into Co. In the process 
of this deformation the particle and hole parameters are crossed and the driving terms 
in the non-Unear integral equations are supplemented by terms depending explicitly on 
these parameters. 

We denote the number of hole parameters by nu and the number of particle parame- 
ters by Up, set n' = nh + Up and define the functions 

which are antiderivatives of the bare energy e(k) and of the kernel K{X) 
define 

1 / l + a„(^|a) \ 
'^^^ = -2^^\ 1 + aoW J- 

A straightening of the contours in (21) leads to 

\n{an{yC))=i%{nh-np) - {yi+N /2-M-nh + np)2r[-'^&{X) 

"h r A.J 

+ £e(X-X))-£e(X-Xp-/ ^K{X-n)\n{l + a„(M)). (79) 

j=l ;=i J^o ^711 

This defines an n'-parametric family of functions, depending on {7jj} and {A.^}- The 
individual functions o„ are then determined by the subsidiary conditions 

o„(?i^) = a„(?i^) = -l, j=l,...,nH, k=l,...,np (80) 

fixing these parameters to a discrete set of values. 
We assume that for sufficiently small temperature 

N/2-M-nh + np = (81) 

and use this condition in our low-temperature analysis without further mentioning. 
Then the eigenvalue ratios take the form 

p„(0|a) = ^«exp{££(X5)-££(Xp-/ dXe(X)z(X)|. (82) 

L;=l j=l J So ) 

The amplitudes (45) consist of a determinant factor times 

Ai°)(a) = expj / ^ !„(',+ 1 . (83) 



. We further 
(78) 
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The determinant factor is treated in appendix D. Inserting (24) into (83) and straighten- 
ing the contours we obtain 



A'n\a)=exJ- f dX [ diuz{X)e'{X-iu)z(ju) 

xexpH/ dXziX) yctHX-X^A-Ycthi^-X) 
I -/eo L.ti Pi 

xexpj-y dXz{X) £(cth(X-X^ + ri)+cth(X-X^-ri)) 

- £ {cth{X -X'j + r\) + cih{X -X']-ri)) 

7=1 



(84) 



Here Cq is a contour infinitesimally close to Co and inside Co- The equation holds for 
= 2M. For ^ 2M additional phase factors may arise. 

8.2 Low-temperature analysis of the non-linear integral equations 

From now on we restrict ourselves to the computation of quantities related to the 
generating function. Then M = N/2. We introduce 



eQ{X) = h- 



47(1 -A2) 
ch(2X)-A 



and 



Mo(A,) = -rin(ao(X + iY/2)) , u{X) = -rin(o„(A, + iY/2|a)) 



(85) 



(86) 



where we suppressed the index n and the dependence on a in the definition of u. 
It follows from (79) that 



u{X)=eo{X) + T 



"p m 
iTiia! + £ e(X - + iY/2) - £ e(X - X) + iY/2) 



7=1 



7=1 

dn 



+ t[ I^KiX-^^)ln(l+.-n, (87) 



'Co -17/2 27li 

where we introduced a' = rja/iTl. A similar equation without the contribution propor- 
tional to r in the driving term holds for uq, 



uo{X) = eoiX) + T f ^K{X-ij)ln(l+e-'^). 



'eo-iY/2 27ti 



(88) 
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Since 6 is bounded on the contour Co — iy, the terms in square brackets in (87) may be 
neglected compared to £o(A,) when T becomes small. Thus, u and mq have the same 
zero temperature limit £. The integrals in (87) and (88) vanish on those parts of the 
contour on which Ree > and are nonzero on their complement. We claim that this 
complement is an interval [—Q,Q\ on the real axis. Indeed, if e satisfies the linear 
integral equation 



where Q is determined by the condition e(±!2) = 0' then e behaves similarly as the 
driving term. It is real and even on M, has at most the two zeros ±2 and is negative on 
[—2, Q\- One can convince oneself that Ree is positive on M — iy which is the lower 
part of the integration contour in (87), (88). This is obvious in the special cases when 
we know e explicitly (e.g. for A = or for /i = 0) and can at least be verified on a 
computer for general values of the parameters A and h. Hence, e satisfies (87) and (88) 
up to corrections at most of the order of T. The function e is the well-known dressed 
energy function. 

Our argument becomes rigorous by employing the following 'generalized Som- 
merfeld lemma' which also allows us to obtain the first and second order temperature 
corrections below. 

Lemma 1. Let let u,f be holomorphic in an open set containing a contour Cj,, and let 
fbe bounded on C„. Let v = Rcm, w = ImM. Assume that v has exactly two zeros Q± 

on C„ separating C„ into a part between Q and on which v is negative and a 
remainder on which v is positive. Assume that zip G Z such that w{Q±) = 2npT. 
Assume that is oriented in such a way that Q- comes before Q+. Then (for T > 0) 



We present a proof of this lemma in appendix C. The lemma is designed in such a 
way that it can be directly applied to the integral equations (87), (88). For small T the 
functions u and mq are close to e and therefore satisfy the requirements of Lemma 1. 
Hence, 




(89) 




(90) 



u{X)=eo{X) + Tri{X) + 



inT^\K{X-Q+) K{X-Q_) 



12 [ u'{Q+) u'iQ.) 




(91) 



where u{X) = u(k) — InipT and 




n (X) = 27ti(a' - + £ e(X - + iY/2) - £ e(X - X) 



+ iY/2). 



(92) 



An equation similar to (91) holds for uq. 
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If we neglect the (D{T'^) terms, equation (91) is already a linear integral equation 
for u. In the next step we resolve the T-dependence of Q±. We make the ansatz 

u{X)=e{'k) + Tui{X) + T\2i'k) + 0{T^), Q± = ±Q + Tdl^ + 0{T^). (93) 
Using that u{Q±) = InipT and e(±!2) = we immediately obtain 



Q± = ±Qt 



-T + 0{T'-), 



(94) 



where we agree upon ui{X) = ui (k) — 2nip. 

Inserting (94) into (91) and expanding up to the second order in T we obtain 



u{X) = eo(X) + Tn{X) + T^r2{X)+ ^K{X-iu)u{n) + 0{T^) , (95) 

J-Q 2ni 



where 



171 



4£'(G) 



Then the hnearity of (95) imphes that wi and U2 satisfy 



U2(k) = r2{%) + ^K{X-p)u2{p). 



(96) 

(97a) 
(97b) 



We can express these functions by means of a set of standard functions that appear 
in the description of the ground state properties of the model. We introduce the dressed 
charge function Z, the root density function p, the dressed phase (|) and the resolvent R 
as solutions of the linear integral equations 



/"2 du 



m 

^(X,v) 



e(X + iY/2) [Q dp 



2ni 
2ni 



+ 



+ I ^^Kik-p)i^{p,v), 



/. 

dp 

-Q 



-Q 2ni 
dp 



K{X-p)p{p), 



R{X,v) = KiX-v) + l_^ ^K{X-p)R{p,v) . 
Comparing (97) and (98) we conclude that 

Ml (k) = 2jti 



(a' - p)Z{X) + 1^<^{X,X']- iY/2) -l^<^{XXj- iY/2) 



U2{X) 



171 



71^ 



71^ 



(98a) 
(98b) 
(98c) 
(98d) 

(99a) 
(99b) 



The root density p will reappear later when we calculate the low-temperature asymp- 
totics of the correlation lengths. Recall that ln(a„(X + iY/2|a)) = —e{'k)/T — ui(k) — 
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U2{'k)T + 0{T^). Thus, we have obtained the solutions of (79) forM = N/2 and np = nn 
up to corrections of order T^. 

It remains to determine the hole and particle parameters X!j and by means of the 
subsidiary condition (80). In terms of the function m the latter reads 

u{X'} - iy/2) = 2mT{hj - 1/2) , u(kl - iy/2) = 2niT{pk - 1/2) , (100) 

where G Z are such that the TJj are located below Co while the are located 
above Cq. Assuming that the quantum numbers hj, pk of the particles and holes are all 
different and that the X'j and X^ are uniquely determined by (100) we may write 

?.) = iY/2+x,., XP,=iy/2 + yp,. (101) 

We restrict ourselves to excitations close to the Fermi edges ±Q for which 

4. = ±Q + 0{T), y% = ±Q + 0{T) (102) 

and denote the number of holes and particles close to the left and right Fermi edge by 
n^, Hp and n^,n^, respectively. Since we assume that nh = np = ri we have 

"ft -"p ="p -"ft (103) 

Using (102), (103) and the identity (^(X, Q) - -Q) = Z{X) - 1, following from 
(98a), (98c), we infer that 

ui{X) = u^P{X) + 0{T), u[^\X)=2ni{{a'-e-p)Z{X)+£). (104) 

Hence, 

u{xp = T2niT{hf-l/2) = ±e'{Q){x^^TQ)+uf\Q)T + 0{T^), (105a) 

u(y%) = ±2rar(p± - 1/2) = ±s'(Q)(y^^ t Q) + ^ {Q)T + 0(7^) , (105b) 

where the signs in the first equations have to be understood as part of the definitions of 
and p^. It follows that 

yt. = Q^ + ^)ipt-y^)+^iT')- (106b) 
Our definition of particles and holes imphes that 

lmx±<0, lm);±>0. (107) 

Equation (106) together with the constraint (107) determines the roots of the 
function 1 + an{X\a) close to the Fermi points to first order in T. To this order they are 
equidistantly located on two parallel lines intersecting the canonical contour So at a 
right angle in iY/2 ^ Q. These roots can be occupied by particles or holes. The picture 
is very much the same as for free Fermions. The parameter £, introduced in (103), 
is the excess of particles over holes at the right Fermi edge which, in conventional 
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language, is describing the 'number of Umklapp processes' involved in the excitation. 
The analogy with free Fermions (the case A = 0) also implies that p is the number of 
roots below 2+ if 2+ is above Co (minus the number of roots between 2+ and So if 
Q+ is below). For p / we cannot naively place particles or holes between Q± and 
Co, since this would be inconsistent with our derivation. For 7^ we have to take into 
account that we might cross particle or hole parameters in the course of deforming the 
contour Co into the contour C„ required by Lemma 1 . In this case the function r\ in 
(92) has to be modified accordingly. In the following we shall avoid such complication 
by simply assuming that p = 0. We shall see that this is already enough to derive the 
leading large-distance asymptotics of the longitudinal correlation functions for small 
temperatures and shall come back to the general case in a forthcoming publication. 

Since p = by hypothesis, we suppress the bar in u, ui etc. in the following except 
in some general statements that are valid for any p eZ. Setting p = means that the 
points Q± +i'Y/2 are located between the roots closest to Co, whence llmtif ((2)| < 7t. 
Using (104) with p = and setting Z = Z{Q) this condition is equivalent to 

Rea' 1 „ Rea' 1 

<^<-— ^ — 7 + — ^- (108) 



Z-i-1 2(1 -Z) Z-i-1 2(1 -Z) 

The latter inequality restricts the range of i for given a', Z. It is easy to see that 
1 /V2 < Z, < 1 for < A < 1. Thus, (108) is always satisfied for £ = -1,0, 1 as long 
as Rea' is small enough. As we shall see, these values determine the leading large- 
distance asymptotics of the correlation functions. Obviously the inequality (108) can 
be satisfied by any £ G Z if Z, is close enough to 1 or if we choose a' appropriately. For 
this reason we keep £ G Z general in our subsequent analysis. 

For the low-temperature analysis of the eigenvalue ratios in the next subsection 
we need u\ up to the first order in T. We insert (104), (106) into (99a) and use 
dy<^(k,v) = R(k,v) /2ni to obtain 

2jtir 



+«(x.-e) -E('.74)-|("74)]}+o(^^)- 009) 

Equations (99b) and (109) together with the corresponding linear integral equations 
then determine ln(o„(X|a)) up to the order T. Using these equations we can readily 
calculate the leading order low-temperature expansion of the eigenvalue ratios. 

8.3 Low-temperature analysis of the eigenvalue ratios 

We start from the representation (82) of the eigenvalue ratio. In order to expand the 
sums in the exponent we use (106). Then 

£ E{x';-) - X Eix';) = -2m{Q+iy/2) 



27iire((2 + iY/2) 



h^-iihj + Pj-l) 



+ 0(7^). (110) 
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It remains to calculate the integral 



/ u(a-) 

/ dXe(X)z(X)= / g;e(X + iY/2)In ^+"' ^ 
■/eo yeo-iY/2 2jii 



(111) 



For this purpose one can employ Lemma 1 as well as the formulae of the previous 
section. Introducing the Fermi momentum 



and the sound velocity 



we obtain altogether 



rQ 

kp =11 d^p(?i), 
J-Q 



2np{Q) 



p„(0|a) =^"exp|-2i(a'-f)ytF 



2nT 
vo 



(a' -efZ^-e^+Y, {hj + Pj-i) 
j=i 



(112) 



(113) 



+ 0{T^). (114) 



8.4 Low-temperature analysis of the amplitudes 

The low-temperature analysis of the amplitudes Al^ '' (a), eqn. (84), is the most tedious 
part of the calculation. Those terms in the expression for the amplitudes that depend 
explicitly on the hole and particle parameters can be expanded using (106). We obtain 
for the respective leading terms 



n(^^e-(^))^^(a.e-(^))^^ 



(-1)^ 



n -2 



(115a) 



n M^'j - ^i) -K)\[ n - K) 

^ sh(X)-Xf+Ti)sh(X)-Xf-Ti) ^ |^ sh(^^+26)sh(^^-26) ^| 111^^^ 



sh^(Tl) 



The integrals in the exponent in (84) which have a non-singular integrand can be 
evaluated by means of Lemma 1 which together with (94) implies the following 



Corollary 1. For f holomorphic in a finite strip around Co and z defined in (78) 

'•Q dX 

/Co J-Q 



[ dXmziX) = - ^f{X + iy/2)uf\X) + 0iT). (116) 

Jen J-0 2Jtl 
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This can be used to conclude that to leading order in T 

- [ dXz{X)\Y,{cth{x-xp +r\) +cth{x-'k'; -^)) 

- ^ {cth{X -X'} +-(])+ cth{X -X'} --(])) 
+ f dX [ dfiz{X)ctW{X-M + y\)z{M) 



Vsh(T| + 2e)sh(Tl-22)y ^ ' J-Q 1-a '^sh2(X-;U + Tl) 

The remaining integrals involve singularities close to the contour of integration and 
are harder to estimate. Fortunately these integrals are similar to integrals that appeared 
previously in the calculation of the low-temperature asymptotics of the generating 
function of the density-density correlators of the Bose gas in [37]. The following lemma 
is slightly stronger than the corresponding statement in case of the Bose gas in appendix 
B of [37] and can be obtained in the same spirit. 

Lemma 2. Let u and Sj, be subject to the same assumptions as in Lemma L Let X+ be 
located above S„ and X- below. Then the Cauchy-type integral 

Iu{'X±) = f dXcih{X-X±)\n[\+Q-'^^ (118) 

admits a low-temperature expansion whose form depends on the distance of X± from 
the zeros Q± of the real part of u. 

For 6 > define V± = {z G C| |m(z)| < b/2^z close to <2±}- If^± uniformly 
away from Q±, then Lemma 1 applies, and 

Iu{X±) = - dXcih{X-X±f-^ + 0{T). (119) 

JQ- T 

IfX±e V+, then 
IfX±eV-, then 

First of all this can be used to calculate the periodic Cauchy transform 

Lejz](v)= / dpz{p)cih{^-^) (122) 
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when V is close to iY/2 ± Q. Assuming this and setting s = sign (im (k)) we find 



+ ■ 



2ni 



In 



sh(X-Q) 



^4Vln(±j27tir) 
2711 



s\ni 



[r{l/2±suo{X)/2niT)J^''^^- ^ 



This formula is useful for the evaluation of the remaining integrals. Using also (105a) 
and (106) we obtain 



2 f dXz{X) £ cth(X - - L cth(A. - X']) 



2iu 



2nT J 



n 

.i=i 



+ 0{T\n{T)). 

(124) 



Finally the most dehcate term is the remaining double integral with singular inte- 
grand. For this integral we can use again (123) and proceed analogously to appendix B 
of [37]. Then to leading order in T 



+: 



(125) 



where G is the Barnes G-function and C\ is a functional defined by 



Ci[v] 



1 rQ 



dX / d/i 



Q ^ v'{X)vi^)-v{X)v'{M) 



(126) 



th(?i-;u) 

When we insert all the terms (115), (117), (124) and (125) into (84) and also use 
(104) they neatly combine. The amplitude A^^ decomposes into two factors aI"^^ 
pertaining to the left and right Fermi edge: 

Ai°) (i7ta/Ti) = aI"^ (a)Ai+^ (a) , (127) 
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where 



(e)sh(2e) 



) 



71 / 



„e^,e{{h)},{p)}\e[{a-i)Z + i]) , (128) 



for e = ±1, and where we have introduced the shorthand notations 



CM 

and 



2{h{X-ij) 



(129) 



^nun2{{hj},{Pj}\v) 



Ul<j<k<ni {hj - hkfX\\<j<k<n2{Pj-Pkf 



(130) 



Thus, in the low-temperature limit we have have obtained an entirely explicit description 
of the 'universal contributions' aI"^ to the amplitudes in terms of the dressed charge 
and the dressed energy. Let us recall that the above formulae (128)-(130) are valid for a 
certain class of excitations characterized by (i) «p = nt, and (ii)x%j% = ±Q + 0{T). 

The Fredholm determinant part of the amplitudes is a fraction with two Fredholm 
determinants in the numerator and two Fredholm determinants in the denominator. In 
order to calculate the zero temperature limit of the determinants in the denominator 
we recall that the weight functions (l + ' iXj) ' and (l + a„ ^ (X|a)) ^ play the role 
of the Fermi functions for our system. In the zero temperature limit they turn into the 
characteristic functions of the 'interval' i'^/2 + [—Q,Q\ on Cq. Thus, 

lim det {l-3C| = lim det{l-X|= det {\-%\. (131) 

The expression on the right hand side depends only on the anisotropy parameter and on 
the magnetic field. For the proof we note the estimate [41] 



|det(l-A)-det(l-S)| < ||A-fi||iexp(||A||i + + 1) , 



(132) 



where ||A||i = tr|A| and which is valid for trace class operators A, B [42]. 

The zero temperature limit of the determinants in the numerator is more technical. 
In order to obtain it we need to assume certain general properties of the low-temperature 
behaviour of zeroes of certain functions (see appendix D). Setting 

^m = ^. . _^o.u......n. (133) 
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I _ g±27ti(a'-£)Z(X) 



and denoting by r[— a contour which encircles the interval [— Q, Q\ counterclock- 
wise we conjecture that 



lim det |l 



det {l 

ciM£,r[-(2,G] 



(134) 
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With (131), (134) the zero temperature limit of the determinant part of the amplitude is 

^ ^^^dM«,r[-Q,e] { 1 - ^-«} dgtdjgg r[-Q,Q] { 1 - ^a} 

The important point about this expression is that it depends only on £ but not on the 
'quantum numbers' pj, hk of the particle hole excitations. 

In appendix D we justify the above conjecture. We also provide alternative expres- 
sions for the ratios of Fredholm determinants that allow us to compare our formulae 
with those obtained for zero temperature in [25] and to extract the a dependence from 
the determinants. Moreover, we provide expressions in which ratios of determinants 
are replaced by determinants of a different type of operators. 



8.5 Summation of the leading terms 

As we already mentioned above it is possible to sum up the form factors of particle-hole 
type considered in the previous subsections. Summation is possible for each individual 
value of £ by means of the same summation formula as employed in [24] in the context 
of the so-called critical form factors pertaining to the eigenstates of the ordinary transfer 
matrix. Adapted to our notation it reads 

„.nt>0 Pi< - <P„„/!, <-</i„. V ^ / 



np,nh>0 Pi< -<Pnphi<-<hn^ 



G^{\+£-v) e 



•'0 



^ ' I 1 Vn 1 



1-e "0 j 



(136) 



Summing the contributions to the form factors series of the form considered in the 
previous subsections amounts to summing over all sets of mutually distinct quantum 
numbers {^j"}, {p^ } and {/ij}, {p^} under the constraint — n'l^ =n^ —n~ =£ and 
then summing over the allowed values of £ G Z. Using (136) we end up with 

/ ttT / \ 2(^-«)^2.2 
/g2maSH \ _ y ^nam+2ikAe-a)m^(^^^J^(^^ _ ^^ ( /^O \ ^^^^^ 

^ l Vsh(m7tr/vo) ) 

where 

A{x) = ^ ^ . (138) 

(2jip(G)sh(2e))'"^ 

For every finite temperature (137) is a sum over terms which decay exponentially 
with the distance m. For T — > 0+ the nature of the asymptotics changes from exponential 
to algebraic decay, (jir/vo) / sh(mnT /vq) — > 1/m, as is expected at a critical point. Up 
to an adjustment of conventions and notation (137) turns into the sum obtained in [24]. 
In particular, the amplitude '!){£) A{£ — a) is the same as in the zero temperature case. 

Looking at it the other way round, the lowest order effect of switching on the 
temperature is a deformation, 1/m ^ {nT /vq)/ sh{mnT /vq), exactly as expected from 
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conformal field theory by mapping the complex plane to a cyUnder of finite circum- 
ference [10]. Here we have obtained this result directly from an expansion of the 
generating function in terms of form factors of the quantum transfer matrix. 

The series (137) is neither an asymptotic series for small T nor for large m. In [37] it 
was characterized as a sum gathering the 'leading orders for each oscillatory term.' This 
means that for each order of exponential decay we have neglected algebraic corrections 
in T which would contribute lower order terms than the next-order exponentials. In 
addition, we have neglected higher order correction to the inverse correlation lengths 
that would amount to contributions of the form exp (0(mr2)). Thus, (137) provides 
only the leading order asymptotics. On the other hand, each oscillating term appears 
with its coefficient of leading order in T . 

Recall that we have assumed in our above analysis that p = 0. This assumption 
restricts the allowed range of summation over ^ for fixed a to the values determined 
by (108). For any AG (0, 1) and any value of the magnetic field this range includes 
£ = — 1 , 0, 1 if a is small enough. It becomes large for magnetic fields close to saturation 
and for anisotropics close to zero. We beheve that an extension of our analysis io p^O 
will show that the sum may be extended to all I ^"L. This would also restore the 
periodicity in a of the generating function. We shall come back to this in future work. 



8.6 Large-distance asymptotics of the longitudinal correlation functions 

The large-distance asymptotics of the longitudinal correlation functions can be obtained 
from (137) by means of (26), (27). From (138) and from appendix D we know that 

D(0) = l + 0(a2), D(^) = 0(a2)for^7^0, (139a) 

yi(0) = l, a) =>[(€) + 0(a) for ^7^0. (139b) 

Then (27) implies 

{o\) = \-^-^ (140) 

which is the known relation between magnetization and Fermi momentum. Using this 
formula and (26) we obtain to leading order (i.e. by taking into account only ^=—1,0,1 
in (137)) 

2 / /„„ \ 2 



Here we have introduced the shorthand notation D"{\) = d\T){\)\a=Q and used that 
D"(-l) = 2)"(l) and^(-l)=yi(l). 

Depending on whether 2. > 1 or 2. < 1 either the first or the second term determines 
the large-distance asymptotics of the correlation function. It is easy to see from the 
integral equation (98a) for the dressed charge function that 2.>lfor— 1<A<0 and 
2,<lforO<A<l provided that Q>0. Thus, as is well known [1, 39], the low- 
temperature correlation functions are asymptotically always negative if — 1 < A < 0, 
while they exhibit Ikp oscillations if < A < 1 . At the free Fermion point A = 
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both types of behaviour coexist. The amplitude in the first term is proportional to the 
square of the dressed charge and is numerically easily accessible. It was obtained using 
a perturbative scheme around A = in [1]. The amplitude in the second term was 
apparently not known until recently [22]. It is not yet clear how well the Fredholm 
determinant part can be computed numerically. This question deserves further attention. 

Note that we have derived (141) assuming that < A < 1. It remains valid for 
— 1 < A < 0, since the correlation functions in this regime of the anisotropy can be 
obtained from those for < A < 1 by a simple transformation [56] consisting of a sign 
reversal of the coupling / and a conjugation of all local operators on even sites by a^. 

9 The finite-length case 

The description of the spectral properties of integrable finite-length systems can be 
based on similar auxiliary functions as introduced above in order to deal with the 
distributions of Bethe roots of the quantum transfer matrix. In fact it is exactly this 
context in which auxiliary functions and non-linear integral equations first became 
relevant for integrable quantum systems [32]. As we have pointed out elsewhere [4, 12] 
the form of many formulae for correlation functions is the same for the infinite chain 
at finite temperature and for the finite chain in the ground state with no magnetic 
field applied. Here we claim that this also holds for the formulae for the form factors 
and amplitudes obtained above. Moreover, if one is willing to accept slightly more 
compUcated contours, then the finite-length case with non-zero magnetic field can 
be treated in a similar manner. The resulting formulae may be problematic for a 
numerical analysis, but the large-L asymptotics may be carried out in analogy with the 
low-temperature analysis of the previous section. This might be seen as an alternative 
to the approach pursued in [23,25]. 

In the finite-length case we work with the monodromy matrix 

= q"''Ra,L{K^/2) . . ./?,,_L+i(X,ri/2) (142) 

acting on the tensor product of an auxiliary space with index a and 2L quantum spaces. 
Here 

K = i(^/ri (143) 

and (j) is a real flux parameter. 

The corresponding transfer matrix = Tra{Ta(}C)} generates the Hamiltonian of 
the XXZ chain, 

H = U,h{y\)r\x\/2)t\^/2)=J f (a5_ia^ + a^_ia^ + A(a}_ia5.-l)), 

7=-L+l 

(144) 

with twisted periodic boundary conditions, 

a^_i = ^--^^ai^'^^S b = x,y,z. (145) 

Since TaQC) satisfies the Yang-Baxter algebra relations (6), the eigenstates and 
eigenvalues of tQC) can again be constructed by means of the algebraic Bethe ansatz. 

The only essential difference is in the pseudo vacuum, which in this case is |0) = (q) 
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and, consequently, in the vacuum expectation values a{X), d{X) of the diagonal elements 
of the monodromy matrix which are now given by 

a{X) = q^, d{X) = ^-V^(?i,ri/2) . (146) 

Inserting these vacuum expectation values into the Bethe equations (16) we obtain 
rather different patterns of roots in the complex plane. Still, the definitions (19) and (20) 
remain meaningful, and a non-linear integral equation similar to (21) can be derived for 
the auxiUary functions. 



ln(o„(?i)) = -2(K-M)ri+2iLp(?i)- / ^KiX-^^)ln(l + an{^i)) . (147) 

Je„ 2m 

Here 



is the bare momentum, whose appearance seems rather naturally in view of the fact that 
we have interchanged 'space and time direction' in the underlying six- vertex model. 
With the auxiliary functions (147) the formula (24) for the ratio of two transfer matrix 
eigenvalues remains valid in the finite-length case if we replace A'^ by 2L. 

The general formula (12) for the correlation functions remains valid as well. Due 
our conventions for the monodromy matrix (142) the operators and eigenvalues on the 
right hand side have to be evaluated at ri/2 instead of 0, 

_ (yo|Tr{0(i)r(Ti/2)}...Tr{QWr(Ti/2)}|>Po) 

The corresponding form factor expansions of the generating function of the longitudinal 
correlation functions and of the transversal two-point function follow as 

Nm-I Nm 

<^''^('")>^= I A„(a)p-(n|a), Ka+^i)^ = i:A„ + (n)p-(n|a) . (150) 

Here the number of overturned spins M must be L for the generating function and L — 1 
for the transversal correlation function. Not only the structure of these asymptotic 
expansions is the same as in the finite-temperature case, but equations (48) and (65) 
for the amplitudes remain valid as well. Note, however, that in the derivation of the 
integral formula (20) for the ratio of the eigenvalues as well as in the derivation of the 
formula for the amplitudes in the asymptotic expansion of the transversal correlation 
functions we have assumed that the spectral parameter is located inside the integration 
contour. This must be kept in mind when sending it to ri/2. In order to deal with a finite 
magnetic field one has to switch to a sector of non-zero eigenvalue of the z-component 
of the total spin. Aq must be replaced by the largest transfer matrix eigenvalue in this 
sector, and the contours have to be changed accordingly. 



10 Conclusions 



In this work we have initiated a study of the finite-temperature correlation functions of 
the XXZ chain by means of the form factors of the quantum transfer matrix. Products of 
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two form factors combine into the amplitudes in the long-distance asymptotics expan- 
sions of the two-point correlation functions. These amplitudes seem to be more natural 
than the form factors themselves. There is no issue about their proper normalization, 
whereas the normalization we chose e.g. for and F+{t,) in (53) seems to be 

somewhat arbitrary. 

For the amplitudes we were able to perform the Trotter limit analytically. The 
resulting formulae (48) and (65) for the longitudinal and transversal cases have a 
remarkably similar structure. They consist of a 'Fredholm determinant part', a 'universal 
amplitude part' , which seems to carry the critical behaviour in the zero temperature 
limit, and a 'factorizing part' by which we mean the product of a-functions in (48) or 
the product of G-functions in (65), respectively. We believe that such type of structure 
persists for more general form factors such as (*F"|C(^i)C(^2)|*i'o)- This conjecture 
is supported by [2, 8, 19], where factorization of correlation functions was proved in a 
rather general context. 

We have demonstrated that our formulae can be efficiently analyzed at least for 
low and high temperatures. This analysis will be continued in future work. We hope 
that the formulae for the amplitudes can also be evaluated numerically. At fixed finite 
temperature a few (typically one or two) amplitudes and the corresponding correlation 
lengths fix the leading asymptotics. Due to the efficiency of asymptotic expansions a 
combination with the exact short-distance results of [3,44,54] would give us numerical 
access to the full correlation functions and would allow us to calculate their Fourier 
transforms which appear in many applications in condensed matter physics. 

An interesting feature of our Fredholm determinants is that they truncate exactly 
for every finite Trotter number (every finite length in the finite-length case). For finite 
Trotter number they are represented by finitely many integrals (see (41)). The Fredholm 
determinants as well as the other factors in the expressions for the amplitudes are 
parameterized by the auxiliary functions o„ whose low-temperature behaviour is rather 
well known. In particular it can be also analyzed for vanishing magnetic field [31]. 
We hope that this fact will help us to analyze the low-temperature asymptotics of the 
full thermal correlation functions for zero magnetic field which is one of the many 
interesting open problems in our approach. 

A further interesting issue to be studied in future work are finite-temperature 
singularities of the amplitudes which may occur e.g. at the cross-over temperatures [13] 
where the asymptotic behaviour of the correlation functions changes from decaying and 
oscillating to monotonically decaying. Finally, we would like to mention the possibility 
to study dynamical correlation functions at finite temperature by means of the form 
factors of the quantum ttansfer mattix. Results for zero temperature, based on the form 
factor approach, were recently obtained in [26]. 
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Appendix A: Amplitudes of the longitudinal correlation functions 

In this appendix we sketch the derivation of the expression (38) for the amplitudes A„(a) 
in the form factor expansion of the generating function of the longitudinal two-point 
functions defined in (30). 

Recall that for the longitudinal correlation functions only excited states of the 
twisted transfer matrix with M = N /2 Bethe roots contribute to the form factor expan- 
sion. We denote these Bethe roots by juj, j = l,...,M, while the Bethe roots of the 
dominant state will be denoted Xj, j = 1 , . . . ,M. Using the scalar product formula (36) 
we obtain 



M 



-7=1 



n 



where we have introduced the function 



detM{e(Xy - iJk) - e{iJk - '>^j)ao{nk)] 



(A.1) 



^ ysh(x-x,-) 

Employing now the Cauchy determinant formula written in the form 

res(|)(Xi:) 



(A.2) 



and the identity 



n '^rf^,=^-h^, ^ 



(A.3) 



(A.4) 



which follows from the Bethe ansatz equations for the /i^, we can simplify the right 
hand side of (A.l). We find that 



sh{fij-fik) 

.i<j<k<MMh-'^k). 

■ M 



■;=i 



A^{Q{Xj-Hk)-Q{Hk-h)^i^k)} 



= det 



shOu -;i^) Jm^^^ "'^ \nk-r[)Q(ki-Hk) + q''^ \nk + V[)Q{nk-h)] 
= deti5^p„'(^j|a)a;(//;|a 



,1+a 



-1-a 



(A.5) 



The contour C in the second equation is chosen in such a way that all Bethe roots 7ij, 
j = 1, . . . , M, are included, but i^j,i^k and /lyt ± t] are excluded. The integral can be 
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calculated, since the integrand is an i7t-periodic fimction of X which decreases to zero 
as ReX ^ ±00.** In the last equation we have inserted the identity 

f (a/;) {q-''r\^^j - Tl) - qy-'i^j+n)) = p-\^Jj\a)aMa) (A.6) 

which follows from Baxter's TQ-equation. Comparing (A.l) and (A.5) and taking into 
account the definition (39) of the kernel Xa we obtain 



detM{5;[-^r(^0C(/i,-/i/,)} 

(A.7) 

A very similar calculation with the roles of Xj and i^k interchanged yields 



1-7=1 



.^q^^d{Xj)e^ 



-7=1 



('I'ol'I'o) 



detM{6i-e^0Ca(X,-X,)} 



(A.8) 

Then the equation (38) for the amplitudes follows by multiplying (A.7) and (A.8). 



Appendix B: Amplitudes of the transversal correlation functions 

In this appendix we derive the expressions (55) and (61) for the form factors and 
F_(^) defined in (52). 



B.l The form factor F+ (^) 

Now the excited state is parameterized by M — 1 Bethe roots still denoted /uj. 
Again a function (|), now defined by 



(B.l) 



will be a useful tool in our calculation. We further use a degenerate variant of the 
Cauchy determinant formula obtained by sending the real part of one of the variables 
(jum in our case) to —00 and picking up the asymptotically leading term, 



n M^^j-^k) 

l<j<k<M-l 

n Mh-h) 

l<j<k<M 



M 

YlTes(\>{Xj) 

1-7=1 



sh(Xi-//i) 



det 

M 



sh(Xi-/iM-i) 
1 



sh(X,M-i"i) 

Sh(XiK-y"Af-ll 
1 



■ (B.2) 



Inserting the determinant formula (36) and the expression (17) for the eigenvalue 
into the definition (52) of F_ (^) and using the above formulae we obtain 



(B.3) 



*A very similar calculation, inspired by [16], was performed in the appendix of [4]. 
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where 



where 



-^7 _ 



, j = l,...,M-l;£=l,...,M, 



.,M, 



(B.4) 

(B.5a) 
(B.5b) 



and where 



{ukY = -Tes^{7ii){q > \iUk-r\)e{h-Hk) + q''^ \iUk + r\Mnk-h)}, (B.6a) 

k=\,...,M-\; e=\,...,M, 



{hmY = - 



l + ao(^) 



res(\>(Xe){e{Xe-^)-e{^-Xe)ao{^)}, £ = \,...,M. (B.6b) 



The determinant on the right hand side of (B.3) can be simphfied in a similar way 
as shown in (A.5) for one of the form factors of the generating function. The four 
cases Cuk, k = I,. . . ,M — I, Cum, Euk, k = I,. .. ,M — I, and Eum have to be treated 
separately. We obtain 



+ ■ 



(Cum) 



£u^ = (^«-i-^i-«)e-«, 

1 q-'-Mm 



(B.7a) 
(B.7b) 



1 



q^nim 



1 



M^-Hj) l + a„(^|a) sh(^-/i,--Ti) l + a„(^|a) sh(^-/i,-+Ti) ' 

(B.7c) 



V l + a„(^|a) l + a^C^la);' 



(B.7d) 



where a„ = 1 /a„ by definition. 

If we directly continue with these matrices we end up with Fredholm determinants 
with kernel %. Having in mind the more involved case of the form factor F_ (^) 
we prefer to switch at this point to the kernel K. For this purpose we subtract q^' 
times (B.7b) from (B.7a) and times (B.7d) from (B.7c) which does not change the 
determinant in (B.3). Then 

{CukY - e^^£u, = bio!„{nj\a)p-' (/iy|a) - e'^-^'^K,.^{nj-Hk) 

= e^'J-^"'{\ky , (B.8a) 

{CuMy - ^^EUM = e^^-^ [cth(^ - ^,-) - cth(^ - - 11) 



l + a„(^|a) 
e^^-^(vM)^ 



■cth(^-^y + ri) 



(B.8b) 
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where j = 1,...,M— 1. These equations define an (M — 1 ) x (M — 1 ) matrix V with 
matrix elements = (v^)^ and an (M— l)-component column vector vm- 

Using (B.8) in (B.3) and expanding the determinant with respect to its last row we 
see that 

= Ae-^ det (y){£uMe^ - (£uie^S . . . ,£uM-ie^^-0^"'vM} • (B.9) 
Here we have assumed that V is invertible, i.e. that 

is non-zero. We do not have a general proof of this assumption, but we checked it with 
several examples for small M. 

It remains to calculate the terms in curly brackets on the right hand side of (B.9). 
Setting X = Vm we find that the coordinates of this vector satisfy the set of Mnear 
equations 

a'MO')p-\Mj\a)x^ = -cthiMj - ^) + - 

This suggests to define a function 

G+ {X, k) = - cth(X - ^) + Y^^^ cth(X - ^ + 11) 

a-l^ ('?|fY'l 

+ \J^^<^H^-^-^)+l.K,.^{^-^,u)x' (B.12) 

which has the following properties: 

G+Oufc,^) = o!M^)9n\^ik\o.)xl' , (B.13) 

and G+(A,,^) as a function of X. is meromorphic inside the contour e„ with a single 
simple pole at X. = ^ with residue — 1 . It follows that G+ satisfies the linear integral 
equation 

G+(X,^) = -cth(X-^) 

+ ^«-ip„(^|a)cth(X-^-ll)+ / dm«(Ai)^:i_a(X-/i)G+(/i,^), (B.14) 



where 



dXp„(X|a) 
+^ ^ 2ra(l + o„(X|a)) 
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Finally 

M~l 

Eums^ - (£uie^> , . . . ,EuM-i^''-')V-'^yM = Eums^ - £ Euj^'^x^ 

= l-^«-ip„(^|a)-(^«-i-^i-«) / dm1{X)G+{X,^) 

= lim G+(X,^) = G;(^). (B.16) 

Here we have used (B.13) in the second equation. Inserting (B.4), (B.IO) and (B.16) 
into (B.9) we arrive at equation (55) of the main text. 



B.2 The form factor F_ (^) 



Here we can proceed as in the previous section by using in (53) the second equation 
(36) with oo on the right hand side. We obtain a formula similar to (55), but with a 
different determinant in the denominator. Multiplying this expression for F (^) with 
(55) we obtain equation (58) for the amplitude A~+(^). In order to obtain the simpler 
and more sjmimetric formula (64) of the main text, however, we have to perform a more 
cumbersome calculation which we shall sketch below. 

Instead of using (36) directly in (53) we first of all move B(^) through the product 
of C operators in (*Fo| to the left, using the well-known Yang-Baxter algebra relations 
(cf. e.g. [15, 36]). This produces a double sum F_ (^) = E^^' ff,m5 where 



■ d(ke)a(km)c(k£,t,)c{i„ X„ 



■M+\ 

n 



1 



M+1 1 

rr — 



k^£,m 



(B.17) 



and Xm+1 = ^ in the various products above. Here we insert the scalar product formula 
(36) on the right hand side. The result can be written in the following form, 

^^■^ = t Iwm^ .S^(-^) det {^(^4 -ri)[e(Ai; -e(X4 -^j)an{hM)]} 

{-ao{Xf:))sign{m-i){-lY-'" 



sh(Xf:-Xm + r\) 

where 4 e {l,...,M+l}\{i,m},k= 1, . . . ,M- 1, £1 < 4 < • • • < Em-i, and 



, (B.18) 



B 



D 



nf=i diXj)e^J detM{dia'oiXj) - K{Xj -Xk)}' 



sh(>.i-pi) 



. g"res(|) ^(w) 
\ sh(XM-//i) 



>ih{Xi-fJM-i) 



g"res(|) '(uM-i) 

sh(X,M-iUM-l) 



(B.19a) 



(B.19b) 
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The function (|) was defined in equation (B.l) above. At this point we see already that the 
prefactor in B which would be somewhat inconvenient in the Trotter limit will cancel 
the corresponding prefactor in A, equation (B.4). What remains to be done is to sum up 
the ti^m into a single determinant. We shall divide this task into several steps. 
Let us define 



where j = 1, 
Then 



{akY = (^{'kk-Ti\){e(jjj-'kk)-e(kk-i^j)a„{Xk\a)) , 
1 

Skj 



(B.20a) 
(B.20b) 



sh(kj-Xk+r[) ' 

M — 1 in (B.20a) and where the remaining indices run from 1 to M + 1 



M 



^P^^ det(D) £ i-lf+'-'sM+i,e det (ai, . . . ,a,_i,a,+i, . . . ,aM) 



l + ao(^) 



detfDf ^> ),...,d( )). (B.21) 



We denote the matrix in the determinant on the right hand side of this equation by X. 
The entries of this matrix can be simplified the same way as in the previous section, and 



+ ■ 



s\v{Xj -Xk-"^) sh(?i; - + ri) 



+ e^isM+i±- (B.22) 



This expression looks appealing as it is a sum of a matrix and a one-dimensional 
projector. One might be tempted to pull the matrix Xl — e^'SM+\.k out of the determinant. 
This is impossible, however, since the matrix is not invertible. The latter fact can be 
easily seen by replacing SM+i,k on the right hand side of (B.21) by zero. A left 
eigenvector with zero eigenvalue has the components res(|)(Xy ). 

In order to resolve the problem we introduce the kernel K at this point and write 



d^t(X) = dti{b{pnQ^j\o)o!Q{Xj)-Ki+aQ^j - Xk)+SM+l,k} , 



1+a 



We further define 



p„(X^|a)a[,(X;) 



(B.23a) 
(B.23b) 

(B.24) 



and assume henceforth that 1 — A' is invertible. This seems to be true for general M. It 
is easy to prove for M = 1, and for larger M we verified it numerically for a number of 
examples. Setting 

w = (l-^)-'v (B.25) 



we conclude that 



det(X) 

M 



M 

n 



det(l-is:) 



M 



1 + -^M+l j-W^ 
7=1 



(B.26) 
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As before we convert (B.25) into an integral equation. It is equivalent to the linear 
equation 



M 



Defining 



we conclude that 



and 



M 



p„{Xj\a)a'QiXj) 



(B.27) 



(B.28) 



(B.29) 



o{X) = 1+ f dm«(;u)i^i+a(^-A')a(A/) . (B.30) 

The function a{X) is a generalized dressed charge function of the type introduced in [5]. 
It allows us to write 



1 + =a.-y^^dm_(X)^j^^^^^) ' 
dXp-\X\a) 



where 



dm«(X) 



27ti(l + ao(X)) 



We finally argue that 



Ooo = lim G(k)=0. 

ReX->+°o 



(B.31) 



(B.32) 



(B.33) 



This can be seen as follows. We could have introduced a parameter, y say, into the last 
column of the matrix D above and at the same time could have multiplied the whole 
sum (B.21) we are considering by l/y. This trivial modification cannot change our 
result, but we would have a term Goo /y instead of Ooo in the final formula, whence (B.33) 
must hold. Inserting (B.26), (B.31) and (B.33) into (B.21) we obtain 



M 



M+l 



e=i l + ao(^) 

In a similar way 



M „/ 



fJiPn\'kj\a) 



det(l-f) [ dm«(?i) 



sh(X-^ + ri) • 

(B.34) 



I = -^^7^ fn det(l -f ) / d.« (^) 



f=\p„'{Xj\a) 



sh(^-?i + ri) ■ 
(B.35) 

To prepare for the calculation of the remaining double sum with ^, m < M + 1 we 
introduce the M xM matrix 



ai ... a^-i aM+i a^+i ... a^ 



(B.36) 
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Then 

M 



1=1 



det(D) 



l + 0o(^) M 

[det([/„) - {-\)^-"'sm,m det (ai, . . . ,a^_i,a;„+i, . . . ,aM)| (B.37) 

L M M—l J 



which still has to be summed over m. 

For the second term in curly brackets we can proceed as above and obtain 

^P" '^^'"^ det(D) f (-if-'^+^Sm.m det (ai, . . . ,a^-i,a^+i, . . . ,aM) 



l + ao(^) -^A-i' 



det(l-f)/" dm«m^^^. (B.38) 
M ' '7e„ sh(ri) 



For the first term we first simplify DUm and perform the summation afterwards. Thus, 

{DU,„)i = bip„iXj\a)a[,iXj) - e^^"^^ {Ki+aih " ^k) - Sm,k) , for A: 7^ m, (B.39a) 
(£>t/^)4 = e^^-^{p„{Xj\aMXj)yj + ^ , (B.39b) 
where 

p„{Xj\aMXj)f = -p„(^|a) (1 + ao® ) cth(X,- - ^) 

+ q»+' cth{Xj - ^ + T1) + ao(^)^-«-i cth(X,- - ^ - Ti) , (B.40a) 



5 = p„(^|a) (1 + oo(^)) - - oo(^)^-"-^ + e^5«,« , 
and hence 



(B.40b) 



det{DUm) = e^"-^ 

M 



.ypnH^.-|a) 



det(l-7S:) 



M 



z'" + 5w'"+£5„,£det 



^=1 



where we defined 



z = (l-f)-V. 
Once again we convert (B.42) into an integral equation. Setting 



(B.41) 



(B.42) 



G-(X,^) = - cth(X-g + '^"'''P" cth(X-^ + Ti) 



l + ao(^) 



+ 



M 



P,7'«|a)z' 



implies that 



^,_ p,(^|a)(l + ao(^))G-(X,-,^) 
p„(Xy|a)a[,(Xy) 



(B.44) 
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and 



G_(X,^) = -cth(X-^) 

+ ?"""Vnn^|a)cth(X-^-Tl)+/" dm'^iM)K,+aiX-n)G.{n,^). (B.45) 

Jen 

Using (B.29) and (B.44) we can calculate the sum on the right hand side of (B.41) 
and can convert it into an integral, 



,-+Jw'" + ^S'^,,detr ' 



= f dm^iX).^ I d^-(^)e^^-(^'^^;y-^-(^'^^"(^^ 

4 ' 4 -^^^ sh(X-;U-Tl) 

+ TTW ie^ IKoiT + "HW 4 sh(^-x+n) 

- / '^^-^^^d^TT^^-^-®/ dm«(X)eMX), (B.46) 

l + ao(q) Je„ sh(/,-^ + ri) Je„ 

where by definition 

Gt{k)= lim G_(X,^). (B.47) 

ReA->+°° 

Notice that in (B.46) we have used the measure dm" when rewriting the finite sums as 
integrals and that we employed (B.33) to simplify the right hand side of (B.46). Using 
a very similar argument as in the derivation of (B.33) we find that 

Gl(^) / dm«(X)eV>^)=0. (B.48) 

With this we can summarize what we have obtained so far. Combining (B.34), 
(B.35), (B.37), (B.38), (B.41), (B.46) and (B.48) we conclude that 

^ e-^nyj/^-"^(A'i)e^- detM{8^^(^/:i+a(X,-X,)} 

X / d««(X)e^ / ^ G_(X,Wa^)-0-(/^,W^) _ 
JQn hn sh(X-/i-ri) 

This is already a quite remarkable formula. The prefactor nicely combines with the 
corresponding prefactor in our expression (55) for The double integral is 

much reminiscent of the double integral representing the inhomogeneous neighbour 
correlation function (Oj^aj ) derived in [4]. In particular, its structure is such that it can 
be factorized by means of the technique suggested in [4, 6]. 
For this purpose we first note that 
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Then 



= f dm«(^)e^| dm«(/i)e^- 



/ 



sh(X-/7-ri) 



e„ 'sh(X-^-Ti) 



Define 



sh(A — /I — T]) q q 

We observe that r{k,iS) — e^^s(k,iLi) is a symmetric function. Hence, we can replace 
r{X,,fj) by e^^.v(A,,/j) under the double integral. Then, further using the integral equations 
(B.30) and (B.45) as well as the asymptotics (B.33) and the elementary identity 



sh(X-^-il) 

= ^«+ie-^+^(cth(X-^-r|) + 1) -^-«-ie^-^(cth(X-^-'n) - l) , (B.53) 
we obtain 



+ ^x3^ / dm«(^)a(X)(-cth(X-^)+^«+ip„U^I«)cth(X-^-il) 

-l + ^«+ip„-i(^|cc))- (B.54) 



Equation (B.48) states that Gl(^) = or the integral is zero. For M = 1 it is easy to 
see that the integral in (B.48) is non-zero. Hence, Gl(^) = 0. We have verified this 
numerically for more examples and in the following assume it to be generally true. 
Then / is reduced to the second term in (B.54). 

Introducing a function G_ as the solution of the linear integral equation 

G_(X,^) = -cth(X-^) 

+ ^«+ip„-i(^|cc)cth(:^-^-il)+ / dm^_{fi)G.{fiX)KiMl^-X) (B.55) 

and using the dressed function trick"''^ we end up with 



where 



(B.56) 



g!(^)= lim G-(X,t,)- (B-57) 
ReX->+~ 



'f'I'Schematically: fi=g\+fKfi and /i = gi + / ^^/2 implies //igi = ffigi- 
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This completes our derivation of equation (61) of the main text. 
We expect that 

limG^(^) = (B.58) 
a->0 

for / to be finite. For the verification note that a = in (B.53) impUes the identity 

Ki (X) = e-^^ [K_ i{X)-q-^+q]+q-^-q (B .59) 

which allows us to replace the kernel Ki by K-i in (B.23a). Proceeding with the 
calculation the function (for a = 0) along with a modified function a will appear 
instead of Gl and a. For these functions an identity similar to (B.48) can be derived 
which implies (B.58). 

Appendix C: A generalized Sommerfeld lemma 

In this appendix we provide a proof of Lemma 1 . Let F be an antiderivative of /. Using 
integration by parts we obtain 



'J d?i/(?i)ln(l+e-'^) 

=n„(2)(f(e,)-f(2_))+rdX^-«"'« 

JQ- I 



+ e T 



= T dXf{X)ln(^l+e'^^ + J^^ dkF{X)u'(X) 

= T dXf{X)\n(\+t^^- dXf{X){u{X)-2%ipT). (C.l) 
Setting uQC) = u{X) — InipT it follows that 

A/ = r^ dXf{X)ln(^l+e-'^^+ J^' dXf{X)u{X) 

= ry^ dX/(X)ln(l+e- ''t'^' Y (C.2) 

Then the real part of the exponent on the right hand side is always negative except at 
Q±. This means that the leading contribution to the integral for T — > comes from the 
(infinitesimally small) vicinities of these two points. In order to quantify the leading 
contribution we fix 5 > small enough. Since u and / are holomorphic we can deform 
the contour locally in a small vicinity of Q± into contours j\. such that w(?i) = l%pT 
for X. G and v(A,) = ±6 at the boundaries of /£. Note that v is monotonic on j\, 
since it has simple zeros at (2±. It follows that 

AI = T f dXf{X)\n(l+e-^)+0{T°°). (C.3) 
We parameterize Jj.hy x = v{X) <^X = (x). Then 



TV f{Q-) 
6 u'iQ.) 



+ (D{T^). (C.4) 
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Here we have used that v~^(0) = Q- and that ln(l +e~l^l) is even. At we can 
perform a similar calculation with the only difference that v{X) is ascending in the 
direction of the contour, whence the sign will be positive. This completes the proof of 
the lemma. 



Appendix D: Properties of the Fredholm determinants 

In this appendix we obtain alternative forms of the Fredhohn determinants or their ratios 
and discuss the zero-temperature Umit of the determinants appearing in the numerator 
of the determinant part of the amplitudes. For the sake of brevity we restrict ourselves 
to the longitudinal case. 



D.l On the zero-temperature Umit 

In order to calculate the T ^ 0+ limit, we note that the TQ-equation implies (see [5]) 

1 ^ ^-XX + ri) 

p„(?i|a)(l + oo(X)) ~ p„{X\a){q^<^{X-y\)-q-»<^{X+y\)) 

+ m (Di) 

where (|) is the ratio of Q-functions defined in (A.2). The function on the left hand 
side is the weight function in the measure dm". The function (|)(X + ri)/p„(A,|a) has no 
poles inside S„. We shall assume that the contour of the Fredholm determinant with 
measure dm" can be deformed into a contour vi \ containing no zero of q^(^{X — Tl) — 
q^"-<^{X + r\) and no pole of the kernel, without changing the value of the determinant. 
In the Fredholm determinant with this deformed contour we can use (D.l) to replace 
the measure dm" by 

In fact, we expect most of the zeros of ^"(|)(X. — ti) — ^~"(|)(X. + t]) to be located close 
to ibTi, where the essential singularities of the functions (|)(XibTi) appear in the Trotter 
Umit. Thus, vi ^ might be similar to S„. 

The weight function in the new measure is easily re-expressed in terms of the 
auxihary functions, 

m l + a„(X|a) 1 ^^^^ 



qy{X-^)-q-y(X + ^) p„(X|a)(l + ao(X)) l-a„(X|a)/oo(X) ' 

We rewrite the factors on the right hand side in a form that is appropriate for performing 
the zero temperature limit. Using the non-linear integral equation (21) we obtain 



an(^|«) 
ao(k) 

while (24) implies that 

l + a„(^|a) 



:^-2«exp|^ dAi^:(Ai-X)z(Ay)|, (D.4) 



p„(X|a)(l + oo(X)) 



:^-"exp|^ dA/e(Ay-X)zw|, (D.5) 
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where X is now outside C„. 

We now assume that we may shift the contour such that its upper part is 
slightly above So, uniformly for all T. Clearly the poles of the function (]) do not prevent 
us from shifting the contour, but we must further assume that — r]) — ^^"(|)(?i+ri) 
has no zeros close to Co for T — > 0+. For the zero temperature limit of the Fredholm 
determinant with measure dM" it is then enough to perform the hmit in (D.4) and (D.5) 
for X away from Co- First straightening the contours in (D.4) and (D.5) and then using 
corollary 1 as well as (104) and (98a) we obtain 



f dfiK{fi - X)z{n) = 2ra(a' -£)- 2ra(a' - e)Z{X - iy/2) + 0{T), (D.6a) 

[ d^ie(ji-X)z{M) = -{a'-i) di^e(j^ + iy/2-X)Zi^) + 0{T) (D.6b) 
Je„ J-o 



l&n J-Q 

which fixes lim7'^o+ d^- (^) for X G t\^\ 

A similar reasoning works for dm" . The identity 

9n{X\a) q"-^-^{X + V[)pn{X\a) 



\ + an{X\a) q-^^-\X-^)-q^^-^{X + ^) 



A ^ ^ CD 7) 

1 (^ - Ti ) - + n ) 



inspires the definition of the measure 

dX (Sf-\X) 



dM«(X) 



2'K\q-^^-\X-X\)-q^^-^{X + r[) 

d?ip„(?i|a)(l + oo(X)) 



27ti l + o„(X|a) 1 - oo(i^)/a„(X|a) 



(D.8) 



which replaces dm" (X) on an appropriate contour T^,^\ We assume it can be kept away 
from Cq. Then eqs. (D.6) determine the zero temperature limit of dM" on F^^'. 
Our claims can be summarized as follows. Define 

±ma!±(a'-t)!%6ne>(n-X)Z{n) 



Then 



lim ml (X) = dMl (X - iy/2) . (D. 10) 

r->-o+ 



hi order to obtain the formulae in the main text it remains to shrink the contours F^"^' to 
a contour F[— 2, 2] + iY/2, encircling [—Q,Q]+iy/2m counterclockwise manner and 
to shift this down to the real axis. 

Note that Z{X) is real on the real axis. Therefore, in accordance with our as- 
sumptions, the denominator in (D.9) is nonzero on the real axis, as long as a' has a 
non- vanishing imaginary part. 
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D.2 Relation to previous determinant formulae 



Using some of the results of the previous subsection we can rewrite the determinants in 
the numerator and establish a relation with the determinants obtained in [25, 37, 38]. 
For this purpose we first note that by construction 



det{l-OC^a}= det {1-3C^4. 



(D.ll) 



Using eqs. (D.4) and (D.5) and inserting the definition (122) of the periodic Cauchy 
transform, we rewrite dM" as 



dM«(X) = T^ 



27ti e'^Le„[z]{X±r\) _ g-2aQTLe„[z]{XTr\) 



(D.12) 



Due to (41) and (47) we may replace X-a{'>^) by K_aQ^) -q^ + and Xa{^) by 
^a(— ^) — + <?" in the determinants in (D. 11). 
We define P = q-^''. 



^:p(X) = cth(?i-ri) - pcth(X + ri) 



and two new kernels 



1 c^^n[m{KpiX-iu)-K^{Q-iu)) 
27ti e^G" H (^-'l) - pe^e„ W 

1 e-'^e„M(/')(^p(X-^)-^p(X-e)) 
27ii e-^enlzK/i+il) _pe-^e„[z](/'-'n) 



Then clearly 



det |1-X^„}= lhn^_det {l + t/(^)} 



In order to compare with [37, 38] we note that 



lim e-^e„M(e±Ti) ^ ^Ie„dXz{X) ^ ^^'^^ ^ ^-i 

6— >— 0° ReX— >— ~ 



whence 



and therefore 



lim 



Z7±l(l-P) 



det {1+f/f^} 



(D.13) 

(D.14a) 
(D.14b) 

(D.15) 

(D.16) 
(D.17) 

. (D.18) 



Using exactly the same reasoning as in appendix A. 3 of [22] one can show that the 
expression in square brackets on the right hand side is independent of 6. Hence, 



det {1-X^a} = 



r ^\ rufl ^ det {l+f/i^H. (D.19) 



Introducing two independent parameters 9i and 02 for the two cases, the determinants 
take the same form as in [37, 38] (see e.g. eqs. (2.14)-(2.18) in [37]). Taking the limit 
r — > 0+ and setting 61 = —Q, 62 = (2 we reproduce the expressions obtained in [22,25]. 
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D.3 Dependence on the twist parameter 

Equation (D.19) is important, because it allows us to understand the a dependence of 
the determinants. Recalling that |3 = q~^^ we infer that det(j;„a g {1 — X=pa} = 0(a) 
unless n = 0. 

This statement remains true for T — > 0+ unless C, = Q. The case ^ = for T — )■ 0+ 
needs a more careful treatment (compare [22, 38]). Consider (133) and (134) for £ = 0. 
Recall that the contour T[—Q, Q\ encircles the interval [—Q, Q\ counterclockwise. Thus, 
we may replace the contour by [—(2, Q\ if we replace at the same time the weight by 

^"27ti(l-e±27cia'ZW) 1,^ ^ )■ ^^-^"^ 

Here X.+ denotes the Mmit from above and X._ the Mmit from below. This replacement is 
justified since the right hand side of (D.20) can be recast as 

dM«m = —( Sh(^- - Q) sh(^ + g - n) ^ ^"'\±ma'±a' Jg^d^efa-X-)(ZM-Z) 
' 27tiVsh(X_+0sh(A,-e-Ti)y' 

(D.21) 

and since for small a' the singularity of the latter expression can be integrated over. We 
further conclude that Um«->od^±(^) = dX./27ti. 

For the ratio of determinants D{t) defined in (135) eq. (D.21) impUes that 

D(0) = l + 0(a^), (D.22) 

which is needed for calculating the longitudinal correlation functions from the generat- 
ing function. 



D.4 A quotient formula 

This subsection is of auxiliary character. We would like to point out that the ratios 
of Fredhohn determinants in the main text can be combined into single Fredhohn 
determinants. This may turn out to be useful e.g. in the low-temperature analysis of the 
model for vanishing magnetic field. 
Consider, for instance, the ratio 

detM{§^'^3C„(X,-X,)} 

D_ = ^— — ^ (D.23) 

detM{5^-Hjfe3C(X,-A,)} 

occurring in eq. (38). We would hke to write it as a single Fredholm-type determinant 
of the form 

D_ = det-i (ei + wi , . . . , Cm + wm) , (D.24) 

M 

where the e^t denote the canonical unit vectors, = 8^ and where the w^; are to be 
determined. Comparing (D.23) and (D.24) we see that they must satisfy the equations 

w^a[,(?i,-) = %a{^j - Xk)p-\'kk\o:) - OCiXj - Xk) 

M 
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These can be turned into an integral equation. Defining 

F(X,v) = 3Ca(?i-v)p;i(v|a)-X(X-v)+ / dm'^iM)XaQ^-l^)Fii^,v) (D.26) 

we conclude that = F(kj,Xk) / cio(^;)- Thus, 

2)_=det-i|5^ + ;%M| ^ det-i{l+^}- (D.27) 
In a similar way we obtain the alternative representation 

2)_=det|5^+^^%M\ ^ det{l+F}, (D.28) 

where 

F(X,v)=3<:(X-v)-0<:a(X-v)p;;i(v|a)+ /" dm(;u)3C(X-/i)F(/i,v). (D.29) 
Similar expressions can also be obtained for the other ratio of determinants in (38). 
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